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SUMMARY 

A description of the finite element Implementation of Robinson's unified 
viscoplastic model Into the General Purpose Finite Element Program - MARC, Is 
presented. To demonstrate Its application, the Implementation Is applied to 
some uniaxial and multlaxlal problems. A comparison of the results for the 
multlaxlal problem of a thick Internally pressurized cylinder, obtained using 
the finite element Implementation and an analytical solution. Is also pre- 
sented. Excellent agreement obtained confirms the correct finite element 
Implementation of Robinson's model. 


INTRODUCTION 

Modern technological development In aerospace, nuclear, and many other 
Industries has resulted In Increasing use of materials at unprecedentedly 
higher temperatures. At high temperatures there Is an Increased propensity 
for Inelastic deformation mechanisms of creep and plasticity. Realistic and 
rational designs of structural components operating at such elevated tempera- 
tures must rely on accurate constitutive model descriptions of these Inelastic 
deformations In structural materials. 

Conventional constitutive models which treat creep and plastic strains as 
separate noninteracting entitles are Incapable of dealing with observed creep- 
plasticity Interactions. The recent years have witnessed, therefore, a con- 
certed effort by numerous researchers to develop constitutive models where the 
Inelastic strain Is not separated artificially Into time-independent plastic 
and time-dependent creep components. These models - called the unified visco- 
plastic models treat all Inelastic strain (plasticity, creep, relaxation, 
etc.) as a unified and time dependent quantity and thus, automatically Include 
any Interactions. A good general overview of numerous viscoplastic models 
developed In recent years and their predictive capabilities Is given by Walker 
(ref. 1) and Llndholm et al. (ref. 2). 

Unified viscoplastic models express the Inelastic strain rates as a func- 
tion of the current values of stress, temperature, and certain Internal state 
variables. The models also Include evolution (or growth) laws, whereby the 
rates of Internal state variables are connected to the current values of these 
variables, stress, and temperature. The viscoplastic models proposed by dif- 
ferent Investigators employ different functional relations. Krleg (ref. 3) 
has, however, shown that most of these models possess the same basic skeletal 
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framework which generally Includes two state variables - one scalar and the 
other tensorlal to account for Isotropic and kinematic hardening, 
respectively. 

In order to make these models more realistic, an attempt Is made to 
Incorporate as much material science as possible Into the continuum model. As 
a result the mathematical structure of these models Is, In general, very com- 
plex. Not only Is the system of constitutive equations associated with these 
unified models highly nonlinear but also mathematically "stiff." Analytical 
tools commonly used for the solution of classical plasticity or creep problems 
are therefore rendered Inapplicable for use with unified models. Many practi- 
cal problems - such as those arising In the analysis of hot gas-path components 
of gas turbine engines and rocket engines for reusable space propul slon, sys- 
tems - Involve complex geometries as well as complex thermomechanical loading 
histories. For application of viscoplastic models to such Intricate problems 
It Is necessary that these models be Implemented In General Purpose Finite 
Element Codes such as MARC (ref. 4). 

During the last several years the NASA Lewis Research Center has mounted 
an extensive In-house and external effort to Implement some of the more com- 
monly used viscoplastic models Into the Finite Element Code - MARC. The visco- 
plastic models due to Walker (ref. 1) and Bodner (ref. 5) have already been 
Implemented Into the MARC and results of some nonlinear structural analyses 
using these Implementations are available In references 1 and 6. As a contin- 
ued effort In this direction, the task of Implementing Robinson's model 
(ref. 7) Into MARC was undertaken. This report provides a brief description 
of this Implementation and also presents the solutions of some Illustrative 
problems . 


ROBINSON'S MODEL 


Robinson's model (ref. 7) Is based on the concept of flow potential. The 
flow and growth laws for the Internal variables are derived from this flow 
potential. The material behavior Is elastic for all the stress states within 
the f low potential , and It Is viscoplastic for all the stress states outside 
the flow potential. Assuming the displacements (and strains) to be small, the 


total strain rate 
(Including creep, 
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The material Is assumed to be Isotropic and Hooke's law Is used to relate the 
elastic strain rate c ^ j to the stress rate, ctij l.e. 


• el 1 + v • u • 

c 1j = E" °1j ~ E °kk 6 1j ’ 

where E Is the Young's modulus, v Is the Poisson's ratio, and 6 ^ j Is 
the Kronecker delta. Following the Cartesian tensor notations, the repeated 
subscripts In equation (2) (and elsewhere) Imply summation. 
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The devlatorlc stress, S^j, is defined as: 
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The growth law governing the evolution of Internal state variable, o^j, 
(which accounts for the kinematic hardening), Is given by: 

®1J = h < tt kl > G 1 j " r ( a kl ) a 1 j • (4) 


The form of this law Is based on the well known Bal ley-Orowan theory 
(ref. 8), which states that the high temperature deformation of materials takes 
place under the Influence of two competing mechanisms denoted by the two terms 
In equation (4). The first term denotes the hardening process with accumula- 
ted deformation and the second term, a recovery or softening process proceed- 
ing with time. Under the steady-state conditions, these two mechanisms balance 

each other and consequently, aij = 0. 

The flow law governing the Inelastic strain rate, e^, Is written as: 
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In equation (5) •Dj denotes the effective stress, which Is defined as: 

Elj = 5 U - a 1j • ( 6 > 

Particular forms for the functions f(F), h(a|<i), and r(a|<-|) are taken as: 
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The parameters y, H, n, m, 0, and G 0 are temperature Independent 
constants. Effects of temperature are Incorporated In the model through the 
temperature-dependent constants y and R using the following relations: 

y = y exp (-©i) , ’’ 


and 


R = 9.0xl0 -8 exp ( 62 ) 


( 12 ) 


In which 

0 1 = (23.8 0 - 2635) ' e) ’ (13) 

and 

e 2 ■ 40 000 (aTT - s) • < H > 

Here, e denotes the absolute temperature In Kelvin (K). The scalar 
state variable * In the above equations denotes the scalar threshold 

(Bingham) stress. Inelastic strain rate, , vanishes for all the values of 

second Invariant J 2 below k. The state variable * accounts for the 
Isotropic hardening (or softening). Following Robinson (ref. 7), and for the 
material considered herein (2-1/4 Cr-1 Mo steel), k has been taken as a 
temperature-independent constant. But, depending on the material, equations 
for the evolution of the Isotropic state variable and for Its dependence on 
temperature, may be necessary. An evolutionary equation for k may be found 
In reference 9. The Inequalities In equations (5) to (9) define boundaries 
across which the flow law and growth law change form dlscontlnuously . 


DESCRIPTION OF MARC - PROGRAM 

Constitutive models are most advantageous If they can be used on a prac- 
tical scale for the solution of problems faced by Industry. The viscoplastic 
model developed by Robinson was therefore Implemented Into the General Purpose 
Finite Element Program - MARC (ref. 4). The MARC - program employs 
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sophisticated Integration algorithms and advanced finite element formulations, 
and Is specially tailored to fit the requirements of a nonlinear structural 
analysis. 

Robinson's model was Implemented Into MARC program by Incorporating all 
the material nonlinearity Into an Initial load vector and treating It as a 
psuedo body force In finite element equilibrium equations. Since the consti- 
tutive equations associated with the viscoplastic model are, In general, 
"stiff," the subincrement technique Is used to form the Incremental constitu- 
tive equations corresponding to the finite load Increment. 

The subincrement technique essentially consists of splitting the finite 
load Increment Into a number of equal subincrements and Integrating the con- 
stitutive equations of viscoplastic models over small subincrements to obtain 
an accurate representation of the Incremental constitutive equations over the 
finite load Increment. The explicit Euler forward difference method (with a 
self-adaptive Integration strategy) was employed to Integrate the constitutive 
equations over the subincrements. The technique Is found to work efficiently 
and accurately, even for large finite element load Increments, provided the 
subincrements are small enough to preserve the stability of Euler forward 
method. However, It Is difficult for the user to select efficient sublncre- 
mental steps, and there Is a considerable Incentive to use as few subincrements 
as possible without affecting the stability of constitutive differential equa- 
tions associated with the model. 

To familiarize the reader with Its operation, a brief Introduction to 
MARC - Program Is appropriate. This Introduction, taken from references 1 and 
10, Is presented In the following paragraphs. 

"The principal of virtual work may be used to generate the MARC nonlinear 
equilibrium equations governing the Incremental response of the structure to 
an Increment of load. In evaluating the nonlinear structural response of a 
component, the program assumes that the load history Is divided Into Incre- 
mentally applied loads. Each load step Is sequentially analyzed as a linear 
matrix problem using an appropriate stiffness matrix and load vector. Although 
each load step uses linear matrix methods to solve the Incremental equations, 
the Incremental equations themselves are nonlinear since the load vector will 
depend on the displacement Increment obtained In the solution of Incremental 
equilibrium equations. 

The principle of virtual work may be written, for applied external point 
loads , or displacement u^. In the form: 

«*{«i dV = 6u{p i , (15) 

where the Integral extends over the volume, V, of each finite element and the 
summation sign extends to all the elements In the structure. 


In equation (15) the virtual displacement vector 6 U ^ Is related to the 
virtual strain vector through the relationship 


6e 1 " B 1j 6 uj or Ae 1 ‘ 6 uj B 1j 


(16) 
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where Is the strain displacement matrix and the superscript T denotes 

transposition. Since A^l Is an arbitrary virtual displacement vector, 
equations (lb) and (lb) may be written In the form: 


l/ v b Vj dv * p , • < 

This relation expresses the equilibrium of structure when the applied 
load vector Is P-^ and the stress vector Is a\. If an Incremental load 
AP^ Is applied to the structure and the stress vector changes to + Ao^ , 
the relation expressing the equilibrium of the structure at the end of the 
Incremental load application may be written as: 



Adj ) dV = P 1 ♦ AP 1 


(18) 


Hence, the relation expressing the equilibrium of the structure during the 
application of the Incremental load vector AP^ Is obtained from equations (17) 
and (18) by subtraction In the form: 



Adj dV 


AP j 


(19) 


Ihe MARC code allows the user to Implement very general constitutive rela- 
tionships Into the program by means of the user subroutine HYPELA. Within this 
subroutine, the user must specify the values of the elasticity matrix 
and the Inelastic stress Increment vector A£i In the Incremental vector 
constitutive relationship: 

A^ = D (Ac A jCi A0) - Af Jj . (20) 

The Inelastic stress Increment vector Afi Is computed In HYPELA using the 
constitutive relationships of the viscoplastic model. 


In equation (20), a denotes the coefficient of thermal expansion and 
Is the vector Kronecker delta symbol. 
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j 


1 If 0 < j < 3 

0 If 3 < j < 6 


( 21 ) 


For the viscoplastic models, "the Incremental Inelastic stress vector A<^ 
depends In a highly nonlinear manner on the Incremental strain vector Ae^. 
Since Aei = A u j , the Incremental Inelastic stress vector Af 1 depends 
In a highly nonlinear manner on the nodal displacement vector Auj , so that 
AC 1 = AC( A u j ) • 


Substitution of equation (2) Into equation (19) produces the Incremental 
equilibrium equations for MARC In the form: 




AP , 


AR 1 + 
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ac j 


dV * 
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A0 dV 


( 22 ) 
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where K^j Is the elemental elastic stiffness matrix defined by the relation: 

■ / v B !k°k« 8 « dv • < 23 > 

The vector AR^ Is the residual load correction vector or out-of-equl 1 Ibrlum 
force vector from the preceding load Increment: 

4R i ■ p i - LJ v 6 ij°j dV • ‘ 24 > 

which Is added to the current Increment In order to restore the structure to 
equilibrium. The nonlinearity In the Incremental equilibrium relationship, 
defined In equation (22), arises because the Inelastic stress Increment vector 
ACj depends nonllnearly on the displacement Increment vector A U j . Values of 

and Afj appropriate to the current Incremental load step are returned 
to the main program by subroutine HYPELA and the Incremental equilibrium rela- 
tions In equation (22) are solved by successive Iterations. 

The solution of the Incremental equilibrium equations In (22) Is accom- 
plished within the MARC code by the following algorithm. At the start of the 
Increment the user subroutine HYPELA Is entered to determine the elasticity 
matrix D^j and the Incremental Inelastic stress vector AC). On entry to the 
subroutine, the Input consists of the strain Increment vector Ac}, the tem- 
perature Increment A0, the time Increment At over which the Incremental 
external load vector AP^ Is applied to the structure, and the values of 
stress, strain, temperature and viscoplastic state variables at the beginning 
of the Increment. Since the Incremental strain vector, Ac^ = A U j, can only 
be accurately determined after the solution to the Incremental equilibrium 
relationship In equation (22) has yielded the correct Incremental solution vec- 
tor A u j, the strain Increment vector Ac^ Initially used to generate the 
Inelastic stress Increment vector AC) must be estimated. The Initial esti- 
mate for Ac] Is assumed to be the value obtained for Ae^ In the preced- 
ing Increment. On exit from subroutine HYPELA the elasticity matrix and 

the estimated Inelastic stress Increment vector AC) are passed Into the main 
program. After the values of and AC) are obtained for each Integra- 

tion point In the structure, the Incremental equilibrium relationship In 
equation (22) Is assembled and solved for the Incremental node displacement 
vector a u j . The Incremental strain vector, <Sc^ = B^j A u j , Is then computed 
and compared with the Initial guess for A used to generate the Inelastic 
Incremental stress vector Afj. If this Incremental strain vector Is equal, 
within a user specified tolerance, to the Incremental strain vector used to 
compute Afj In the assembly phase, the solution Is assumed to have converged. 
Otherwise, the updated strain Increment vector, obtained from the solution of 
the equilibrium relations In equation (22), Is passed Into subroutine HYPELA, 
a new vector, ACj , Is computed and the equilibrium equations resolved to 
yield an Improved value of A u i and Ae^. The process Is repeated until 
the value of vector Ac) on the assembly phase Is equal, within a user 
specified tolerance, to the value of the vector Ae^ on the solution phase. 
After convergence Is achieved, the temperature, stress vector, strain vector 
and viscoplastic state variables are updated by adding the Incremental values 
generated during the current Increment to the values of these variables at the 
beginning of the Increment. The program then passes on to the next load 
Increment where the process Is repeated.' 1 
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1MPLEMEN1A1 ION OF ROBINSON'S MODEL IN MARC 


The MARC - Program provides the user a convenient way of Implementing and 
Integrating the constitutive equations of Robinson's model through the user 
subroutine HYPELA. The details of subroutine HYPELA may be found In ref- 
erence 4. The subroutine HYPELA written for Robinson's model and the self- 
adaptive Integration scheme used In It, follow essentially the same structure 
as given In reference 10 for Walker's model. 

The discontinuous boundaries In Robinson's model should be smoothed to 
facilitate numerical computations using the FEM. The smoothing Is achieved by 
defining a spline function P(x) on the Interval (-1,1) as: 


P(x) = 


(1 » x)' 

2 


1 - 

i ; 
0 ; 




-1 < x < 0 
2 ; 0 < x < 1 


2 

x > 1 
x < -1 


(25) 


The function F In equation (1) Is then replaced by a function !F defined as: 


& = P 


j <F> 


(26) 


where the angular brackets denote: 


<x> = 


X 

0 


x > 0 
x < 0 


(27) 


and the weighting function W^ Is to be selected by the user. The use of 

equation (26) results In smoothing the discontinuity In across the sur- 
face, = 0, and also rids the flow law, equation (5), of Inequalities. 


The discontinuities In c^j 
G by G] , 1 .e. 


are removed by first replacing the function 


G > 2G 


o 


4G 


+ G 


G < 2G 


o 


(28) 


and then defining 


G - ( G ] - G 




( 29 ) 
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I quatlons (8) and (9) may now be rewritten as: 


h 


= G 0 


and 


r 



( 30 ) 


(31) 


With these values of h and r, the growth law, equation (6), reduces to a 
single expression with a smooth transition across the surface, S-tja^j = 0, In 
the space and now contains no Inequalities. The weighting function W 2 Is again 
again to be selected by the user. The computations and results presented In 
the following sections used: 

W-, = W 2 = lO- 2 ( ks 1 2 ) = 0.475 (MPa 2 ) . 


NUMERICAL VALUES OF THE CONSTANTS 

The numerical values of the parameters appearing In equations (1) to (14) 
for the alloy 2-1/4 Cr-1 Mo steel were taken from Robinson and Swlndeman, 

(ref. 7). The values are: 

v = 3.61xl0 7 (hr) 

n =4 
m = 7.73 

B = 1.5 (32) 

H = 1.37x10-4 (ksl/hr) 

G 0 = 0.14 

k = 0.82 ( ks 1 ) 

The values of Young's modulus, E, and the Poisson's ratio, v are assumed 
to be temperature dependent. The following polynomial expressions In temper- 
ature taken from Sartory (ref. 11), are used to estimate the values of E 
and v : 

E = 31100.0 - 13.59 T + 0.2505xl0‘ 5 T 2 - 0.2007xl0- 13 T 3 ( ks 1 ) , (33) 

and 

0 = 0.524 + 0.154X10- 3 T - 0.126xl0- 6 T 3 . (34) 

Here, T denotes temperature In degree Farenhelt. 
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APPLICATION TO PROBLEMS 


Uniaxial Problems 

Some Illustrative problems were solved to demonstrate the Implementation 
of Robinson's model Into MARC through the subroutine HYPELA. The values of 
constants listed In the preceding sections were utilized for these computa- 
tions. The thermal and mechanical loadings used to generate the hysteresis 
loops using MARC are shown In figure 1. 

Figures 2 through 4 show the hysteresis loops generated by MARC using 
Robinson's model. The loops shown In these figures are at strain rates of 
0.0004, 0.004, and 0.04/mln, respectively. The strain range Is +0.2 percent 
approximately. The cyclic mechanical loading used In generating the hystere- 
sis loops Is shown In figure 1(a) and the thermal loading Is constant at 
800 °F . These figures exhibit the results for two and one-quarter of a cycle. 
Figure 5 summarizes the stabllzed hysteretlc loops for the three strain rate 
conditions . 

The results for combined In-phase thermomechanical loading, shown In 
figure 1(b), are depicted In figures 6 to 8. These figures show the hystere- 
sis loops for the same strain-rates and strain range as mentioned above. 

Figure 9 summarizes the stable hysteretlc loops for this type of thermomechani- 
cal loading. 

As a further example, the combined out-phase thermal and mechanical load- 
ing shown In figure 1(c), was used to generate the hysteresis loops at strain 
rates of 0.0004, 0.004, and 0.04/mln, respectively. The strain range, as 
before, Is + 0.2 percent. Figures 10 to 12 show the results of these calcula- 
tions for two and a one-quarter of loading cycles. The stable hysteretlc 
loops for the three strain rates are plotted In figure 13. The loops of 
figures 9 and 13 are Identical If the signs of stresses and strains are simply 
reversed. Figures 2 to 13 depict the effects of different strain rates and 
different types of thermomechanical loadings. 

Figure 14 compares the Isothermal hysteretlc loops (1000 °F) generated by 
MARC code with those obtained experimentally by Robinson and Swlndeman 
(ref. 7). The strain rates are 0.0004 and 0.04/mln, respectively and the 
strain range +0.32 percent, approximately. The excellent agreement between 
the theoretical and experimental loops verifies the correct finite element 
Implementation of Robinson's model. 

Multlaxlal Problem 

In order to have a further verification and application of the finite 
element Implementation of Robinson's constitutive equations, the Implementa- 
tion was applied to a multlaxial problem. The problem selected was that of a 
thick-walled cylinder subjected to Internal pressure under Isothermal (800 °F) 
conditions. The Inner and outer radii of the cylinder were taken to be 
0.16 In. (4.06 mm) and 0.25 In. (6.35 mm) and the internal pressure as 3.65 k s 1 
(25.17 MPa). The results of these computations have been plotted In 
figures 15 and 16. 
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Figure 15 shows the hoop (circumferential) stress distribution across the 
wall of cylinder using the Robinson's model and the FE Code - MARC. The curves 
plotted at different values of time show the redistribution of stress as pre- 
dicted by Robinson's model. The hoop strain distribution obtained using 
Robinson's model and MARC has been plotted In figure 16 for different values 
of time. 

In order to verify the results from MARC and to ensure the correct 
Implementation of Robinson's model In the code, the analytical solution for 
the problem of a thick, Internally pressurized cylinder was obtained. The 
following section presents briefly the details of the analytical solution. 


ANALYTICAL SOLUTION FOR THICK-CYLINDER 


In the following, the subscripts r, q>, and z refer to the radial, hoop 
(circumferential) and axial directions respectively. Symbol, a denotes the 
stress and the symbol e, strain. 

Following Saada (ref. 12), the equation of equilibrium may be written as: 


3o <3 a 
r _ <p - r 

3r r 


(35) 


Assuming the total strain rate to be composed of elastic and Inelastic 
components, the constitutive relations In the rate form may be written as: 


* 1 T* /• • \1 »1n 

e r - E [ tf r - * * °z/J + c r 

• if* /* • \1 *1n 

<p E L <p \ z r/J <p 

• i r* /• • yi • i n 

£ z = E [°z - v \°r * + c z 


(36) 

(37) 


(38) 


•In *1n *1n 

In equations (36) to (38), e r , , and c z denote the Inelastic strain rates 

In the respective directions. These rates are obtained using the constitutive 
equations of Robinson's model. The dot over a symbol Indicates differentia- 
tion with respect to time, t. 


The straln-rate-compatlblllty relation Is: 


3c 

<£ 

dr 



r 


o 


(39) 
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and, the boundary conditions of the problem are: 


a - -p at the inner radius r = a , 

r 

c 3=0 at the outer radius r = b . /(in 

r ( 40 

Using equations (35) to (40), and after some algebraic manipulations, one 
may obtain the following expression for the stress rates: 



These equations together with the constitutive equations of Robinson's 
model were Integrated using the explicit Euler forward method with an automatic 
time-step Integration strategy proposed by Arya et al . (ref. 13). The Integra- 
tion yields the stresses and strains at different times and radii of the cylin- 
der. The spatial Integrals In equations (41) to (43) were evaluated using the 
Trapezoidal rule. 

The values of the hoop (circumferential) stress across the wall of the 
cylinder at different times have been plotted In figure 17. Figure 18 exhibits 
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the values of hoop strain at different times In the cylinder obtained using the 
analytical solution. The values of stress and strain from the MARC program 
have also been plotted In these figures for comparison of results. The excel- 
lent agreement between the MARC and the analytical solutions verifies the 
correct Implementation of Robinson's model Into FE Code - MARC. The finite 
element Implementation can therefore be safely employed to carry out the non- 
linear finite element analysis of structures with more complex geometries - 
such as gas turbine airfoils subjected to complex thermal and mechanical 
loadings In aerospace propulsion engines. The results of such analyses will 
be reported subsequently. 


SUMMARY AND CONCLUSIONS 

1. Robinson's unified viscoplastic model has been successfully Implemen- 
ted Into the Finite Element Code - MARC. 

2. Some uniaxial and multlaxlal problems have been solved using the Imple- 
mentation. The results obtained for a multlaxlal thick cylinder problem using 
the Implementation compare very well with the corresponding results from the 
analytical solution and verify the correct finite element Implementation of 
Robinson's model . 

3. The Implementation can now be used for the analysis of more complex 
problems, both In geometry and loadings, from aerospace and other Industries, 
where analytical methods become Inapplicable, or are too cumbersome to apply. 
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(A) ISOTHERMAL LOADING. 




(B) NON ISOTHERMAL (IN-PHASE) LOADINGS. 




(C) NON ISOTHERMAL (OUT-PHASE) LOADINGS. 



FIGURE 1. - CYCLIC THERMAL AND MECHANICAL LOADINGS USED TO FIGURE 2. - HYSTERESIS LOOP USING ROBINSON'S MODEL AT 800 °F 

GENERATE HYSTERESIS LOOPS BY FE CODE-MARC. (STRAIN RATE = 0.0004/min). 
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FIGURE 3. - HYSTERESIS LOOP USING ROBINSON'S MODEL AT 800 °F 
(STRAIN RATE = 0.004 min). 


FIGURE 4. - HYSTERESIS LOOP USING ROBINSON'S MODEL AT 800 °F 
(STRAIN RATE = 0.04/min). 



-.0020 -.0010 0 .0010 .0020 
STRAIN 


FIGURE 5. - HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES USING 
ROBINSON'S MODEL-ISOTHERMAL CASE (800 °F). 



STRAIN 


FIGURE 6. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NONISO- 
THERMAL CASE (STRAIN RATE = 0. 0004/MI N), IN-PHASE THERMAL 
AND MECHANICAL LOADINGS. 


15 





STRESS/ KSI 




FIGURE 7. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NON ISO- 
THERMAL CASE (STRAIN RATE = O.OCWmin), IN-PHASE THERMAL 
AND MECHANICAL LOADINGS. 


FIGURE 8. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NONISO- 
THERMAL CASE (STRAIN RATE = 0.04/MIN), IN-PHASE THERMAL 
AND MECHANICAL LOADINGS. 



FIGURE 9. - HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES USING 
ROBINSON'S MODEL-NON ISOTHERMAL CASE. IN-PHASE THERMAL AND 
MECHANICAL LOADINGS. 



FIGURE 10. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NONISO- 
THERMAL CASE (STRAIN RATE = 0. 0004/MI N), OUT-PHASE THERMAL 
AND MECHANICAL LOADINGS. 
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STRAIN 

FIGURE 11. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NONISO- 
THERMAL CASE (STRAIN RATE = 0.04/min), OUT-PHASE THERMAL 
AND MECHANICAL LOADINGS. 



STRAIN 

FIGURE 13. - HYSTERESIS LOOPS AT DIFFERENT STRAIN RATES USING 
ROBINSON'S MODEL- NON ISOTHERMAL CASE. OUT-PHASE THERMAL AND 
MECHANICAL LOADINGS. 



STRAIN 

FIGURE 12. - HYSTERESIS LOOP USING ROBINSON'S MODEL-NONISO- 
THERMAL CASE (STRAIN RATE = 0.04/min), OUT-PHASE THERMAL 
AND MECHANICAL LOADINGS. 



FIGURE 14. - COMPARISON OF MARC AND EXPERIMENTAL HYSTERESIS 
LOOPS (1000 °F). 
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FI SURE 15. - STRESS DISTRIBUTION IN AN INTERNALLY PRESSURIZED 
THICK CYLINDER USING ROBINSON'S MODEL (FEM-MARC SOLUTION). 
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FIGURE 16. - STRAIN DISTRIBUTION IN AN INTERNALLY PRESSURIZED 
THICK CYLINDER USING ROBINSON'S MODEL (FEM-MARC SOLUTION). 
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FIGURE 17. - COMPARISON OF MARC AND ANALYTICAL RESULTS FOR 
STRESS DISTRIBUTION IN AN INTERNALLY PRESSURIZED THICK 
CYLINDER USING ROBINSON'S MODEL. 


2-8 tr- 


£a 


2.0 


1.2 


M 


A 

O 

□ 

O 


/A 


TIME - 0.0007 hr (MARC) 

TIME = 0.5500 hr (MARC) 

TIME = 20 hr (MARC) 

TIME - 2000 hr (MARC) 

TIME - 0.0007 HR (ANALYTICAL) 
TIME = 0.5500 hr (ANALYTICAL) 
TIME = 20 hr (ANALYTICAL) 

TIME - 2000 HR (ANALYTICAL) 


‘*A 


^AA 


La 






%A 


<3* 




LA 










.16 .19 -22 .25 

RADIUS, IN. 

FIGURE 18. - COMPARISON OF MARC AND ANALYTICAL RESULTS FOR 
STRAIN DISTRIBUTION IN AN INTERNALLY PRESSURIZED THICK 
CYLINDER USING ROBINSON'S MODEL. 
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